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We present a study on the direct determination of the rf mass on the full set of 
SESAM and TxL QCD vacuum configurations with 2 active flavours of Wilson 
fermions, at j3 = 5.6. We observe a definite dependency of the two-loop correlator 
on the topological charge sector. 



1 Introduction 

Flavour symmetric hadronic states are of particular importance for the un- 
derstanding of nonperturbative aspects of quantum chromodynamics (QCD). 
In fact they are expected to reveal important insight into the topological vac- 
uum structure induced by gJuonic self interactions, as borne out in the famous 
Witten-Veneziano formula El that relates the flavour octet /singlet mass split- 
ting, Mq 0, to the topological susceptibility, Xg: m the large- N c limit of the 
theory: 

Af 2 := M* - Ml = 2N fXq /fl , (1) 

with Nf being the number of active flavours and f v the pion decay constant. 

While glueballs have been the target of many ab initio lattice investiga- 
tions over the past fifteen years □ much less attention has been paid to the 
direct calculation of flavour singlet states with fermionic content such as the 
rf meson. The reason is obvious: the extraction of flavour singlet mesons from 



"Upper (lower) case letters refer to masses in physical (lattice) units. 
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Table 1: Lattice parameters and numbers of independent vacuum field configurations and 
stochastic sources, Ng . Last two columns: t-ranges for the single exponential fits to the 
octet and singlet channel correlators. 
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the lattice is even more challenging than the one of glueball states, due to the 
substantial cost of computing quark (rather than Wilson) loop correlators. In 
fact such Zweig-rule forbidden diagrams (in general misleadingly called 'dis- 
connected diagrams') consist of two quark loops connected via gluon lines only 
and need to be calculated in the momentum zero state, which amounts to the 
computationally very expensive evaluation of the trace of the inverse Dirac 
operator. The two pioneering studies in the field used quenched vacuum con- 
figurations, with Wilson fermions and wall sources Q and staggered fermions 
with stochastic sources □. 

This situation will change with the imminent advent of near-Teraflops 
computers that promise both to provide (a) sufficient sampling rates for full 
QCD vacuum configurations and (b) the analysing power to deal with the 
above mentioned loop-loop correlators. Albeit present day simulations are still 
restricted to the case of two active flavours, Nf = 2, it is of great interest to 
fathom state-of-the art techniques to deal with Zweig rule forbidden diagrams. 
First results on flavour singlet masses were presented recently by CP-PACSB 
and UKQCDB. 

In the present paper we present a study on the rj' mass using the full set of 
QCD configurations generated by SESAM (16 3 x 32, 'small lattice') and T^LLl 
(24 3 x 40, 'large lattice'), both at j3 = 5.6 with standard Wilson action (see 
Table |). 



2 Disconnected diagrams 
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2.1 The problem 

We consider the pseudoscalar flavour singlet operator in a flavour symmetric 
theory 

N f 

= ^2<ii{xh5qi(z) , (2) 

i=l 

with Nf flavours. By the usual Wick contraction it leads to the flavour singlet 
propagator in terms of the inverse Dirac operator, M : 

C v ,(0\x) ~ Ar / tr((M" 1 (0|a;)) t M- :L (0|a;)) 

-7V / 2 tr( 75 M- 1 (0|0)) t tr( 75 Af- 1 (x|x)) , (3) 

which is a sum of fermionic connected and disconnected terms with traces to 
be taken in the spin and colour spaces. In the rest of the paper we shall refer to 
them as 'one-loop' and 'two-loop' contributions, respectively. The momentum 
zero projection 

C, r (t) ee (r/(t)7/(0)) conn - (r/(t)r/(0)) disc (4) 

is expected to decay exponentially, ~ exp(— m^t) and reveal the flavour sin- 
glet mass, 171^1. On an antiperiodic (of length T in time) lattice one should 
encounter the usual cosh behaviour at large values of t and T — t 

G(t) -> exp(-m 7) 't) + exp(-m,,/ (T - 1)) . (5) 

From this parametrization local effective masses m^'(t) can be retrieved by 
solving the implicit equations 

G(t+1) = cxp(-m^(t + 1)) + cxp(-m v ,(T - t - 1)) 
G(t) exp(— m v 't) + cxp(— m v > (T — t)) 

where the l.h.s. ratios are taken from the lattice. 

2.2 Diagonal improved stochastic estimator 

Obviously the momentum zero projection embodied in eq. ||requires the knowl- 
edge of the traces, tr(75M _1 (a;|a:)), on each time slice. But their exact evalu- 
ation is far beyond the scope of present computers. 

One therefore takes resort to a stochastic estimate of the disconnected di- 
agrams by introducing an ensemble of N e nonlocal sources, <j)(x) a , a — 1, . . . N e 
with uncorrelated stochastic entries on the lattice sites x £ colour x spin x i? 4 
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Figure 1: Correlation functions on the small lattice at the lightest sea quark mass, n = 0.1575, 
with local sources. Upper (lower) curves refer to one-loop (two-loop) contributions. 



space such that in the limit N e —> oo the ensemble averages show the following 
limiting behaviours 



(d>(x)) 
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x,y 



(7) 
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In our case, we chose noise for the components of the source vectors. By 
solving the Dirac equation 



M(z,x)£(x)=<t>(z) 



(9) 



on the stochastic sources <j)(z) one retrieves a bias free estimate for the inverse 
Dirac operator in the large iV e limit 



(0(2/)^)) = ^Af- 1 (z, a; )(0(y)t0(z)) - M-\y,x) , 



(10) 



and hence for the one-loop term, tr^M -1 ). But beware of the subleading 
terms when computing components of M _1 as we do here (actually in spin 
space); for they might be affected by the dominant diagonal ones in so far the 
latter are not sufficiently rejected by our approximation to the Kronecker S in 
eq. H So we improve the estimate on nondiagonals, Af~ 1 (?/, x) (x ^ y), by 
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Figure 2: r/ propagators (k = 0.158) with local source and sink on the large lattice. 

subtracting out, from the expression on±he l.h.s. of eq. those leading error 
terms by redefining the estimator to be 13: 

+ M-\z,x)(<j>{y?<t>{z)) -(H) 

Note that in this improvement step, the subleading term on the l.h.s. 
again is gained by SET. Throughout this work we shall apply eq. [ll] in order 
to achieve improved estimates to tr^M -1 ). We mention in passing that an 
alternative route -wipuld be to work in a spin explicite mode by using spin 
projected sources HQ. 

2.3 Quality of the signal 

In Tabic [t] we collected the characteristic parameters of our present analysis. 
We have used five different sea quark masses and two different lattice sizes to 
have some control on finite-size effects. While the number of vacuum configu- 
rations varies from 156 to 195, the number of independent stochastic sources 
has been chosen to be 400 (100) for the small (large) lattices. The psfiudoscalar 
and vector mass ratios quoted refer to our final spectrum analysisE3. The re- 
sults presented in this talk are based on local sources only. The errors quoted 
are statistical and have been obtained by jackknifing. 
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Figure 3: Chiral extrapolation of m 2 , (left) and m^i (right), both linear in the quark mass. 
Masses in lattice units. Full squares: pseudoexperimental value, according to eq. [l6| 



As a result we obtain two-loop correlators, Cdi SC (t)i with statistical qual- 
ity illustrated for our lightest quark mass on the small lattice in Fig. [j]. The 
plot also contains the connected correlator, C 8 (t), corresponding to the nonsin- 
glet pseudoscalar meson, in order to demonstrate the reduced accuracy of the 
(Zweig-rule forbidden) two- loop contribution. At the time separation t ~ 10, 
we are faced with a statistical error of order « 25 % on the latter. 

Given this situation it is practically impossible to establish an effective 
mass plateau in the rf channel, C v >(t), (eq. ||). For further analysis we made 
single exponential fits to (t) and C% (t) over the t-ranges listed in Table [l| 
The quality of our fits is illustrated for the lightest quark mass on the large 
lattice in Fig. g. 

3 Physics results 

3.1 Chiral extrapolations 

Because of— the well-known technical limitations of the hybrid Monte Carlo 
algorithm tJ, the SESAM and T%L configurations correspond to two mass 
degenerate light sea quark flavours (Nf — 2), with the unrenormalized mass 
value (in lattice units) 

m q = 1/2( K - 1 -K- 1 ) . (12) 
From our previous light spectrum analysis 0, we quote the lattice spacing 

a' 1 {ki) = 2.302(64)GeV (13) 

and the critical and light quark n values: 

k c = .158507(44) , K light = .158462(42) . (14) 



G 



From our data we cannot decide whether it is mi, or that follows a linear 
behaviour in the quark mass: the two fits shown in Fig. ^| work equally well, 
with x 2 /d.o.f. ~ 0(1). Note that we make no distinction between sea and 
valence quarks, as we choose the quark masses in the fermion loops to equal 
the sea quark masses (symmetric extrapolation in the sense of ref . 113) . 
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Figure 4: The sensitivity of the disconnected diagram, Rq(t) := C ( ji sc (t)/C8(t), when sub- 
sampling the vacuum configurations according to topological charge, \Q\. 

In our Nf = 2 world, according to eq. [| we would not expect to encounter 
the full effect of Zweig rule fordbidden diagrams, and hence we anticipate to 
underestimate the real world rf mass. From the experimental mass splitting 



A/ 2 = Af 2 - Ml , 



(15) 



we therefore compute, in the spirit of the Witten-Veneziano formula eq. |l], the 
'pseudoexperimental' value in the Nf = 2 world: 



M*,{N f 



2/3M 2 + Ml = (716MeV) 2 



(16) 



This value corresponds in lattice units to the full squares in the two alternative 
chiral extrapolations shown in Fig. 0. To compare: our lattice analysis yields 
at the light quark mass 



Af 2 , = (551(85)MeV) 2 and M v > = 615(53)MeV , 
by use of a linear ansatz in m 2 /(m 9 ) and m^/(m g ), respectively. 
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Figure 5: The effective locai octet (w) and singlet (rj 1 ) masses after Wuppertal type smearing 
on source and sink. 



3. 2 Impact of topology 

We expect two-loop diagrams to be sensitive to the topology of the vacuum. 
A simple check is to look for a dependency of the ratio of disconnected and 
connected correlators, R,Q(t), on the size of the topological charge \Q\, as 
determined in ref.Ej. In Fig. [|, we have plotted, for K sea = .1575 and the small 
lattice, the data with cuts according to \Q\ < 1.5 (topOl) and \Q\ > 1.5 (top24). 
We definitely find a dependency of Rq on \Q\. Note that the disconnected piece 
vanishes in the vacuum sector with small values of |Q|! 



4 Discussion and outlook 

The analysis presented here looks rather promising, but leaves room for further 
investigations: 

Smearing. The rf correlator being a difference of one-loop and two-loop 
terms it is of great importance to attain early asymptotics of the flavour singlet 
correlator. But with local sources the connected (flavour octet) correlator does 
require analysis at t- values beyond 12, in view of appreciable excited state 
contributions. Yet for the flavour singlet situation, we had to analyse in the 
t-range 6 to 9, for lack of statistics. Thus it appears very promising to improve 
the ground state overlap by use of source and sink smearing; for this will help to 
establish mass plateaus in the flavour singlet channel and hence reduce errors. 
Indeed, we have observed mass plateaus appearing at t ~ 4 to 5, as illustrated 



8 



in Fig. |5|. In a forthcoming paper, we shall elaborate on this pointliJ. 

Spectral methods. Another possible direction to go is to construct the 
two-loop correlator from the low lying eigenfunctions of the Dirac operator. 
For illustraton we show in Fig. ^ the result of an attempt to saturate the spec- 
tral representation of the two-loop correlator by the 300 lowest eigenmodes of 
(75M)li3. We find very nice agreement between the results from this eigen- 
mode approach (EVA) from SET. This is very encouraging, since in the deeper 
chiral regime EVA will become even more efficient whereas the performance of 
SET will deteriorate. 
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Figure 6: Saturation of the two-loop correlator with the 300 lowest eigenmodes of the op- 
erator (75 M), at the smallest sea quark mass on the small lattice. Symbols: crosses from 
EVA, horizontal lines from SET. 
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